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Abstract 

There has been major progress in recent years in the development of 
improved discretizations of the QCD action, current operators, etc for 
use in numerical simulations that employ very coarse lattices. These 
lectures review the field theoretic techniques used to design these dis- 
cretizations, techniques for testing and tuning the new formalisms that 
result, and recent simulation results employing these formalisms. 



1 Introduction 

Lattice quantum chromo dynamics is the fundamental theory of low-energy 
strong interactions. In principle, lattice QCD should tell us everything we 
want to know about hadronic spectra and structure, including such phe- 
nomenologically useful things as weak-interaction form factors, decay rates 
and deep-inelastic structure functions. In these lectures^] I discuss a rev- 
olutionary development that makes the techniques of lattice QCD much 
more widely accessible and greatly extends the range of problems that are 
tractable. 

The basic approximation in lattice QCD is the replacement of continuous 
space and time by a discrete grid. The nodes or "sites" of the grid are 
separated by lattice spacing a, and the length of a side of the grid is L: 



1 These lectures were given at the 1996 Schladming Winter School on Perturbative and 
Nonperturbative Aspects of Quantum Field Theory (Schladming, Austria, March 1996) 
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The quark and gluon fields from which the theory is built are specified only 
on the sites of the grid, or on the "links" joining adjacent sites; interpolation 
is used to find the fields between the sites. In this lattice approximation, the 
path integral, from which all quantum mechanical properties of the theory 
can be extracted, becomes an ordinary multidimensional integral where the 
integration variables are the values of the fields at each of the grid sites: 

JvA,...e-f Ldt ^ J J] dA»(*i) ...e-°£ L '. (1) 

Thus the problem of nonperturbative relativistic quantum field theory is 
reduced to one of numerical integration. The integral is over a large number 
of variables and so Monte Carlo methods are generally used in its evalua- 
tion. Note that the path integral uses euclidean time rather than ordinary 
minkowski time, where i cuc i = Itmink! this removes a factor of I from the 
exponent, getting rid of high-frequency oscillations in the integrand that are 
hard to integrate. 

Early enthusiasm for this approach to nonperturbative QCD gradually 
gave way to the sobering realization that very large computers would be 
needed for the numerical integration of the path integral — computers much 
larger than those that existed in the mid 1970's, when lattice QCD was 
invented. Much of a lattice theorist's effort in the first twenty years of lat- 
tice QCD was spent in accumulating computing time on the world's largest 
supercomputers, or in designing and building computers far larger than the 
largest commercially available computers. By the early 1990's, it was widely 
felt that teraflop computers costing tens of million of dollars would be es- 
sential to the final numerical solution of full QCD. 

The revolutionary development I discuss here was the introduction of 
new techniques that allow one to do realistic numerical simulations of QCD 
on ordinary desktop workstations or even personal computers. To under- 
stand this development one must understand the factors that govern the 
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cost of a full QCD simulation. This cost is governed by a formula like 



where the first factor is just the number of lattice sites in the grid, and the 
remaining factors account for the "critical slowing down" of the algorithms 
used in the numerical integration. This formula shows that the single most 
important determinant of cost is the lattice spacing. The cost varies as the 
sixth power of 1/a, suggesting that one ought to keep the lattice spacing as 
large as possible. Until very recently it was thought that lattice spacings as 
small as .05-.lfm would be essential for reliable simulations of QCD. As I 
describe in these lectures, new simulation results based on new techniques 
indicate that spacings as large as .4fm give accurate results. Assuming that 
the cost is proportional to (1/a) 6 , these new simulations using coarse lattices 
should cost 10 3 -10 6 times less than traditional simulations on fine lattices. 

The computational advantage of coarse lattices is enormous and will cer- 
tainly redefine numerical QCD: the simplest calculations can be done on a 
personal computer, while problems of unprecedented difficulty and precision 
can be tackled with large supercomputers. In these lectures I explain why 
we now think coarse lattices can be made to work well. I review the tech- 
niques from quantum field theory needed to redesign lattice QCD for coarse 
lattices. These techniques are based on ideas from renormalization theory 
and effective field theory. I begin in Sect. 2 discussing the field-theoretic im- 
plications of discretizing space and time. Then in Sects. 3, 4 and 5 I discuss 
in detail how to discretize the dynamics of gluons, light quarks and heavy 
quarks, respectively. I also discuss the crucial issue of testing and tuning 
improved discretizations. These sections present a mixture of very old and 
very new results. I do not discuss in any detail the numerical techniques 
used in the simulations; these are described in standard texts JD, 0, [J. As 
these are lectures from a school, I include a number of exercises that il- 
lustrate the concepts developed in the lectures. Finally I summarize the 
current situation and prospects for the future in Sect. 6. 

2 Field Theory on a Lattice 

In this section I discuss the factors that limit the maximum size of the lattice 
spacing. Replacing space-time by a discrete lattice is an approximation. If 
we make the lattice spacing too large, our answers will not be sufficiently 
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accurate. For the purpose of these lectures, I define a "sufficiently accurate" 
simulation to be one that reproduces the low-energy properties of hadrons to 
within a few percent, which is very precise for low-energy strong-interaction 
physics. The issue then is: How large can we make a while keeping errors 
of order a few percent or less? 

A nonzero lattice spacing results in two types of error: the error that 
arises when we replace derivatives in the field equations by finite-difference 
approximations, and the error due to the ultraviolet cutoff imposed by the 
lattice. I now discuss each of these in turn, and then focus on the key role 
played by perturbation theory in correcting for finite-a errors. 



2.1 Approximate Derivatives 

In the lattice approximation, field values are known only at the sites on 
the lattice. Consequently we approximate derivatives in the field equations 
by finite-differences that use only field values at the sites. This is com- 
pletely conventional, and very familiar, numerical analysis. For example, 
the derivative of a field ip evaluated at lattice site Xj is approximated by 

^«A^fe) (3) 

where 

A^(x) ^ ^ + °)-^(*-°). (4) 
2a 

It is easy to analyze the error associated with this approximation. Taylor's 
Theorem implies that 

2a/S. x ip{x) = ip{x + a) — ip(x — a) 

= (e ad * -e~ a9x ) tl){x) (5) 

and therefore 

A^= (d x + ^d* + 0(a 4 yjrp. (6) 

Thus the relative error in A. x ip is of order (a/A) 2 where A is the typical 
length scale in ip{x). 

On coarse lattices we generally need more accurate discretizations than 
this one. These are easily constructed. For example from (||) it is obvious 
that 

^ = A^-^A^ + 0(a 4 ), (7) 
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which is a more accurate discretization. When one wishes to reduce the 
finite-a errors in a simulation, it is usually far more efficient to improve 
the discretization of the derivatives than to reduce the lattice spacing. For 
example, with just the first term in the approximation to d x ip, cutting the 
lattice spacing in half would reduce a 20% error to 5%; but the cost would 
increase by a factor of 2 6 = 64 in a simulation where cost goes like 1/a 6 . On 
the other hand, including the a 2 correction to the derivative, while working 
at the larger lattice spacing, achieves the same reduction in error but with 
a cost increase of only a factor of 2. 

Equation ^ shows the first two terms of a systematic expansion of the 
continuum derivative in powers of a 2 . In principle, higher-order terms can be 
included to obtain greater accuracy, but in practice the first couple of terms 
are sufficiently accurate for most purposes. Simple numerical experiments 
with a 3 -accurate discretizations like this one (see below) show that only 
three or four lattice sites per bump in ip are needed to achieve accuracies 
of a few percent or less. Since ordinary hadrons are approximately 1.8 fm 
in diameter, these experiments suggest that a lattice spacing of .4 fm would 
suffice for simulating these hadrons. However QCD is a quantum theory, 
and, as I discuss in the next section, quantum effects can change everything. 

Exercise: Show that 

0=A$^-£(a«)Vo(« 4 ) (8) 

where 

(2) 

Note that A x and A^ can be used to construct lattice approximations for 
derivatives of any order: for example, 

A(?V = A x AfV> = A^AsV 

= 6gV + 0(o 2 ) (10) 

A<*ty — (Ai 2 )) 2 V 

= 0^ + O(a 2 ) (11) 

Exercise: Use 

^( x ) = e -* 2 /2- 2 (12) 

with (j = 2cz as a sample function to test our discretizations of d x ip an d d 2 ip. 
Evaluate the a 2 errors of the simplest discretization in each case, and the a 4 
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errors of the more accurate discretizations. Show that the a 2 errors are of 
order 10-20%, while the a 4 errors are only a few percent or less. Note that 
A^ 2 ) is more substantially more accurate than A^. 

2.2 Ultraviolet Cutoff 

The shortest wavelength oscillation that can be modeled on a lattice is one 
with wavelength A m i n = 2a; for example, the function ip{x) = +1, —1, +1 . . . 
for x = 0, a, 2a . . . oscillates with this wavelength. Thus gluons and quarks 
with momenta k = 2ir/X larger than ir/a are excluded from the lattice theory 
by the lattice; that is, the lattice functions as an ultraviolet cutoff. In simple 
classical field theories this is often irrelevant: short- wavelength ultraviolet 
modes are either unexcited or decouple from the long-wavelength infrared 
modes of interest. However, in a noisy nonlinear theory, like an interacting 
quantum field theory, ultraviolet modes strongly affect infrared modes by 
renormalizing masses and interactions. Thus we cannot simply discard all 
particles with momenta larger than n/a; we must somehow mimick their 
effects on infrared states. This we can do by modifying our discretized 
lagrangian.0 

To see how we might mimick the effects of k>n/a states on low momen- 
tum states, consider the scattering amplitude T for quark-quark scattering 
in one- loop perturbation theory (Fig. |l]a). The difference between the cor- 
rect amplitude T in the continuum theory and the cut-off amplitude in 
our lattice theory involves internal gluons with momentum k of order n/a 
or larger. This is because the classical theories agree at low momenta, and 
therefore all propagators and vertices agree there as well. Thus the loop 
contributions from low momenta will be the same in T and T^ a \ and cancel 
in the difference. Given that the external quarks have momenta pi that are 
small relative to vr/a, we can expand the difference T — in a Taylor 
series in api to obtain 

T - T (a) = a 2 c(a) u(p 2 )^^u(pi) u(p A )^u(pz) 

+ a 2 c A (a) n7 M 7 5 un7 M 7 5 u 

+ a 4 d(a) (pi — p2) 2 ujfiuuj^u 

+ (13) 

2 The idea of modifying the lagrangian to compensate for a finite UV cutoff is central 
to chiral field theories for pions, nonrelativistic QED/QCD and all other effective field 
theories. The application to lattice field theories was pioneered in the form discussed here 
by Symanzik and is refered to as "Symanzik improvement" of lattice operators H] . 
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Figure 1: One-loop amplitudes contributing to: a) qq^-qq, and b) the quark 
self energy. 

where the couplings c, ca ■ ■ • are dimensionless functions of the cutoff. This 
difference is what is missing from the lattice theory; it is the contribution of 
the k>7r/a states excluded by the lattice.^] The key observation is that we 
can reintroduce this high-fc contribution into the lattice theory by adding 
new interactions to the lattice lagrangian: 

+ | a 2 c A (a) ^7^75^ V ; 7 M 75 ip 
+ a 4 d(a)^-D 2 7 M V?7 M V' 

+ (14) 

These new local interactions give the same contribution to T(qq — > qq) in the 
lattice theory as the k> ir/a states do in the continuum theory. Although 
these correction terms are nonrenormalizable, they do not cause problems in 
our lattice theory because it is cut off at ir/a. On the contrary, they bring our 
lattice theory closer to the continuum without reducing the lattice spacing. 

This simple analysis illustrates a general result of renormalization theory: 
one can mimick the effects of states excluded by a cutoff with extra local 
interactions in the cut-off lagrangian. The correction terms in 5C^ are 
all local — that is, they are polynomials in the fields and in derivatives of 
the fields — since contributions omitted as a consequence of the cutoff can 
be Taylor-expanded in powers of ap{. This is intuitively reasonable since 
these terms correct for contributions from intermediate states that are highly 
virtual and thus quite local in extent; the locality is a consequence of the 
uncertainty principle. The relative importance of the different interactions 

3 More precisely, this contribution is due to states above the cutoff, and to corrections 
needed to fix up the states below but sufficiently near the cutoff that they suffer severe 
lattice distortion. 
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is determined by the number of powers of a in their coefficients, and that 
is determined by the dimensions of the operators involved: each term in 
$£( a ) must have total (energy) dimension four, and so if the operator in 
a particular term has dimension n + 4 then the coefficient must contain a 
factor a n . In principle there are infinitely many correction terms, but we 
need only include interactions with operators of dimension n + 4 or less to 
achieve accuracy through order (api) n . In practice we only need precision to 
a given order n in api, and there are only a finite number of local operators 
with dimension less than or equal to any given n + 4. 

Initially the k>ir/a contributions appear to be bad news. As we argued 
in the previous section, a 2 corrections are necessary for precision when a is 
large. However the quantum nature of our theory implies that there are 
a 2 terms, due to k> ir/a states, which depend in detail on the nature of the 
theory. Thus, for example, when we discretize the derivative in the QCD 
quark action we replace 

Ipd ■ 7V> -> -0A • 7V> + a 2 d(a) V^A (3) • 7^ H (15) 

where now d(a) has a part, —1/6, from numerical analysis (equation (^)) 
plus a contribution that mimicks the k> n/a part of the quark self energy 
(Fig. 0b). The probl em is that the latter contribution is, by its nature, 
completely theory specific; we cannot look it up in a book on numerical 
analysis. We must somehow solve the k > it /a part of the field theory in 
order to calculate it; otherwise we are unable to correct our lagrangian and 
unable to use large lattice spacings. 

The good news, for lattice QCD, is that k > n/a QCD is perturbative 
provided a is small enough (because of asymptotic freedom). Then correc- 
tion coefficients like d(a) can be computed perturbatively, using Feynman 
diagrams, in an expansion in powers of a s (ir/a). Thus, for example, our 
corrected lattice action for quarks becomes 

C (a) = ^(A-7 + m(a))V 
+ d{a) a 2 ^A (3) • 7^ 

+ c{a)a 2 ^^ip^ip + ... (16) 



where 



d{a) = -- + dia s (ir/a) + ■■■ (17) 
6 



c(a) = C2 a 2 (-K/a) + ■ ■ ■ (18) 



8 



are computed to whatever order in perturbation theory is necessary. In this 
way we use perturbation theory to, in effect, fill in the gaps between lattice 
points, allowing us to obtain continuum results without taking the lattice 
spacing to zero. 

2.3 Perturbation Theory and Tadpole Improvement^ 

Improved discretizations and large lattice spacings are old ideas, pioneered 
by Wilson, Symanzik and others ||. However, perturbation theory is es- 
sential; the lattice spacing a must be small enough so that k ~ ir/a QCD 
is perturbative. This was the requirement that drove lattice QCD towards 
very costly simulations with tiny lattice spacings. Traditional perturbation 
theory for lattice QCD begins to fail at distances of order 1/20 to 1/10 fm, 
and therefore lattice spacings must be at least this small before improved ac- 
tions are useful. This seems very odd since phenomenological applications of 
continuum perturbative QCD suggest that perturbation theory works well 
down to energies of order 1 GeV, which corresponds to a lattice spacing 
of 0.6 fm. The breakthrough, in the early 1990's, was the discovery of a triv- 
ial modification of lattice QCD, called "tadpole improvement," that allows 
perturbation theory to work even at distances as large as 1/2 fm []|, [7j. 

One can readily derive Feynman diagram rules for lattice QCD using 
the same techniques as in the continuum, but applied to the lattice la- 
grangian || . The particle propagators and interaction vertices are usually 
complicated functions of the momenta that become identical to their con- 
tinuum analogues in the low- momentum limit. All loop momenta are cut 
off at kn = ±ir/a. 

Testing perturbation theory is also straightforward. One designs short- 
distance quantities that can be computed easily in a simulation (i.e., in 
a Monte Carlo evaluation of the lattice path integral). The Monte Carlo 
gives the exact value which can then be compared with the perturbative 
expansion for the same quantity. An example of such a quantity is the 
expectation value of the Wilson loop operator, 

W(C) = (0|iRe Tr Pe~ l9 & A " da; |0), (19) 

where A is the QCD vector potential, P denotes path ordering, and C is any 
small, closed path or loop on the lattice. W(C) is perturbative for sufficiently 
small loops C. We can test the utility of perturbation theory over any range 

4 This section is based upon work with Paul Mackenzie that is described in 0. 
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of distances by varying the loop size while comparing numerical Monte Carlo 
results for W(C) with perturbation theory. 

Fig. |2] illustrates the highly unsatisfactory state of traditional lattice- 
QCD perturbation theory. There I show the "Creutz ratio" of 2a x 2a, 
2a x a and a x a Wilson loops, 

, fW(2ax2a)W(axa)\ , . 

plotted versus the size 2a of the largest loop. Traditional perturbation theory 
(dotted lines) underestimates the exact result by factors of three or four for 
loops of order 1/2 fm; only when the loops are smaller than 1/20 fm does 
perturbation theory begin to give accurate results. 

The problem with traditional lattice-QCD perturbation theory is that 
the coupling it uses is much too small. The standard practice was to express 
perturbative expansions of short-distance lattice quantities in terms of the 
bare coupling aq at used in the lattice lagrangian. This practice followed 
from the notion that the bare coupling in a cutoff theory is approximately 
equal to the running coupling evaluated at the cutoff scale, here a s (ir/a), 
and therefore that it is the appropriate coupling for quantities dominated 
by momenta near the cutoff. In fact the bare coupling in traditional lattice 
QCD is much smaller than true effective coupling at large lattice spacings: 
for example, 

Qtiat = a v (n/a) - 4.7 ay H (21) 

< \a v {-K/a) forO.lfm (22) 

where cty(q) is a continuum coupling defined by the static-quark potential, 

V m (q) = -4nC F ^. (23) 

Consequently ai at expansions, though formally correct, badly underestimate 
perturbative effects, and converge poorly. 

The anomalously small bare coupling in the traditional lattice theory 
is a symptom of the "tadpole problem". As we discuss later, all gluonic 
operators in lattice QCD are built from the link operator 

rather than from the vector potential An. Thus, for example, the leading 
term in the lagrangian that couples quarks and gluons is ipU^^ip/a. Such 
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Figure 2: The \22 Creutz ratio of Wilson loops versus loop size. Results 

from Monte Carlo simulations (exact), and from tadpole-improved (new) 

and traditional (old) lattice perturbation theory are shown. 
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a term contains the usual ipgA ■ jip vertex, but, in addition, it contains ver- 
tices with any number of additional powers of agA^. These extra vertices 
are irrelevant for classical fields since they are suppressed by powers of the 
lattice spacing. For quantum fields, however, the situation is quite differ- 
ent since pairs of A^'s, if contracted with each other, generate ultraviolet 
divergent factors of 1/a 2 that precisely cancel the extra a's. Consequently 
the contributions generated by the extra vertices are suppressed by powers 
of g 2 (not a), and turn out to be uncomfortably large. These are the tadpole 
contributions. 

The tadpoles result in large renormalizations — often as large as a factor 
of two or three — that spoil naive perturbation theory, and with it our in- 
tuition about the connection between lattice operators and the continuum. 
However tadpole contributions are generically process independent and so 
it is possible to measure their contribution in one quantity and then correct 
for them in all other quantities. 

The simplest way to do this is to cancel them out. The mean value uq 
of | Re Tr consists of only tadpoles and so we can largely cancel the 
tadpole contributions by dividing every link operator by uq. That is, in 
every lattice operator we replace 

U,(x) - ^ (25) 
u 

where uo is computed numerically in a simulation. 

The wo's cancel tadpole contributions, making lattice operators and per- 
turbation theory far more continuum-like in their behavior. Thus, for exam- 
ple, the only change in the standard gluon action when it is tadpole-improved 
is that the new bare coupling axi is enhanced by a factor of 1/uq relative 
to the coupling a\ a t in the unimproved theory: 

cm = (26) 

Since Uq < .6 when a > .lfm, the tadpole- improved coupling is typically 
more than twice as large for coarse lattices. Expressing axi in terms of the 
continuum coupling ay, we find that now our intuition is satisfied: 

«ti = ay^/a) — .5 aty + ■ ■ ■ (27) 
« a v (ir/a). (28) 



Perturbation theory for the Creutz ratio fl20|) converges rapidly to the 
correct answer when it is reexpressed in terms of a-ri- An even better result 
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Figure 3: The critical bare quark mass for Wilson's lattice quark action ver- 
sus lattice spacing. Monte Carlo data points are compared with perturbation 
theories in a theory with no tadpole improvement (T.I), tadpole-improved 
gluon dynamics, and tadpole-improved quark and gluon dynamics. 

is obtained if the expansion is reexpressed as a series in a coupling constant 
defined in terms of a physical quantity, like the static-quark potential, where 
that coupling constant is measured in a simulation. By measuring the cou- 
pling we automatically include any large renormalizations of the coupling 
due to tadpoles. It is important that the scale q* at which the running cou- 
pling constant is evaluated be chosen appropriately for the quantity being 
studied || f?J. When these refinements are added, perturbation theory is 
dramatically improved, and, as illustrated in Fig. ^, is still quite accurate 
for loops as large as 1/2 fm. 

This same conclusion follows from Fig. ^ which shows the value of the 
bare quark mass needed to obtain zero-mass pions using Wilson's lattice 
action for quarks. This quantity diverges linearly as the lattice spacing 
vanishes, and so should be quite perturbative. Here we see dramatic im- 
provements as the tadpoles are removed first from the gluon action, through 
use of an improved coupling, and then also from the quark action. 

The Creutz ratio and the critical quark mass are both very similar to 
the couplings we need to compute for improved lagrangians. Tadpole im- 
provement has been very successful in a wide range of applications. 
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Exercise: Derive the Feynman rules for lattice 4> 4 theory where 



T. ^ 



4! 



(29) 



In particular show that the <f> propagator is 



,2 



G oc 5 (30) 

(2sin(a PAl /2)) 2 + a 2 TO 2 

2.4 An Aside on "Perfect Actions" 

Although not necessary in practice, it is possible to sum all the terms in 
the a 2 expansion of a continuum derivative using Fourier transforms: for 
example, we replace 

dli>(xi) -> -J2 d p{( x i- X j)^( x j) ( 31 ) 

3 

where 

d$(x t - Xj ) = J ^i? e(\p\<n/a). (32) 

We do not do this because S3 (xi — Xj) falls off only as — Xj\ 2 for 

large \xi — xA, resulting in highly nonlocal actions and operators that are 
very costly to simulate, particularly when gauge fields are involved. So it is 
generally far better to truncate the a 2 expansion than to sum to all orders. 
However this analysis suggests a different strategy for discretization that we 
now examine. 

(2) 

The slow fall-off in the all-orders or perfect derivative d pi is caused by 
the abruptness of the lattice cutoff (the 9 funtion in (^)). Wilson noticed 
that d (2) can be made to vanish faster than any power of l/\xi — xA if a 
smoother cutoff is introduced. In his approach the lattice field at node Xi 
is not equal to the continuum field there, but rather equals the continuum 
field averaged over a small volume centered at that node: for example, the 
lattice field at site X\ in one dimension is 5 

i rxi+a/2 

*i = - / ij(x)d X (33) 

a Jxi-a/2 



5 The smearing discussed here is conceptually simple, but better smearings, resulting 
in improved locality, are possible. In || the authors use stochastic smearing in which \f; 
equals the blocked continuum field only on average. This introduces a new parameter that 
can be tuned to optimize the action. 
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30 ^ e i^n( P )v#) (34) 

where, in the transform, smearing function 

s 2simW2) 
pa 

provides a smooth cutoff at large p. The lattice action for the smeared 
fields is designed so that tree-level Green's functions built from these fields 
are exactly equal to the corresponding Green's functions in the continuum 
theory, with the smearing function applied at the ends: for example, 

pXi+a/2 rXj+a/2 

(ViVj) = dx dy (V>(x)V(y)> (36) 

JXi—a/2 Jxj—a/2 

30 jg ffr) (37) 

oo 2-7T p* + m z 

By rewriting the propagator in the form 

* /a dp e i P ( Xi - Xj ) n 2 (p+27rn/a) 
- w /o2tt V {p+2nn/a) 2 + m 2 ' { ' 

we find that the quadratic terms of this lattice action are 

S« = IX^dgfc-^m)*^, (39) 



where dpa? (a;. 



" /a e IK^-%)/\- n 2 (p+27rn/a 



7 7 ? • ( 4 °) 

Va2vr \^ (p+27rn/a) 2 +m 2 y ^ ; 

(2) 

It is easily shown that, remarkably, dp a (x{— Xj, m) vanishes faster than any 
power of l/\xi — Xj\ as the separation increases. So the derivative terms in 
such an action are both exact to all orders in a and quite local — that is, 
the classical action is "perfect." Again, this is possible because the lattice 
field *$>i is obtained by smearing the continuum field, which introduces a 
smooth UV cutoff in momentum space rather than the abrupt cutoff of fl32|) . 

The complication with this approach is that it is difficult to reconstruct 
the continuum fields ip{x) from the smeared lattice fields ^j. Consequently 
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nonlinear interactions, currents, and other operators are generally compli- 
cated functions of the and much more difficult to design than with the 
previous approach, where the lattice field is trivially related to the con- 
tinuum field. This is an important issue in QCD where phenomenological 
studies involve a wide range of currents and operators. It also means that 
in practice the "perfect action" used for an interacting theory is not really 
perfect since the arbitrarily complex functions of *$>i that arise in real perfect 
actions must be simplified for simulations. Although these difficulties grad- 



ually are being overcome for QCD |LC], 11], I will not discuss this interesting 
approach further in these lectures. 



2.5 Summary 

Asymptotic freedom implies that short-distance QCD is simple (perturba- 
tive) while long-distance QCD is difficult (nonperturbative). The lattice 
separates short from long distances, allowing us to exploit this dichotomy to 
create highly efficient algorithms for solving the entire theory: k>n/a QCD 
is included via corrections SC to the lattice lagrangian that are computed 
using perturbation theory; k<n /a QCD is handled nonperturbatively using 
Monte Carlo integration. Thus, while we wish to make the lattice spacing a 
as large as possible, we are constrained by two requirements. First a must be 
sufficiently small that our finite-difference approximations for derivatives in 
the lagrangian and field equations are sufficiently accurate. Second a must 
be sufficiently small that tr/a is a perturbative momentum. Numerical ex- 
periments indicate that both constraints can be satisfied when a ~ 1/2 fm 
or smaller, provided all lattice operators are tadpole improved. 

It is important to remember that most any perturbative analysis in QCD 
is contaminated by nonperturbative effects at some level. This will certainly 
be the case for the couplings in 5C. However we now have extensive expe- 
rience showing that such nonperturbative effects are rarely significant for 
physical quantities at the distances relevant to our discussion. Should a sit- 
uation arise where this is not the case (and one surely will, some day) we may 
have to supplement the perturbative approach outlined in this section with 
nonperturbative techniques. Such situations will not pose a problem if they 
are relatively rare, as now seems likely. Furthermore, as we discuss later, 
simple nonperturbative techniques already exist, should they be needed, for 
tuning the leading correction terms in both the gluon and quark actions. 
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3 Gluon Dynamics on Coarse Lattices 



In this section I discuss first the construction of accurate discretizations of 
the classical theory of gluon dynamics. I then discuss the changes needed 
to make a quantum theory, and illustrate the discussion with Monte Carlo 
results. Finally I discuss actions for anisotropic lattices, and the tuning and 
testing of gluon actions. 



3.1 Classical Gluons 

The continuum action for QCD is 

S = ( d^x^TiF^x) (41) 

where 

F^ u = d^A u - d v A^ + ig[An,A v ] (42) 

is the field tensor, a traceless 3x3 hermitian matrix. The defining charac- 
teristic of the theory is its invariance with respect to gauge transformations 
where 

F^^n(x)F^n(x)l (43) 

and Q(x) is an arbitrary x-dependent SU3 matrix. 

The standard discretization of this theory seems perverse at first sight. 
Rather than specifying the gauge field by the values of A^(x) at the sites of 
the lattice, the field is specified by variables on the links joining the sites. 
In the classical theory, the "link variable" on the link joining a site at x to 
one at x + afi is determined by the line integral of A^ along the link: 



U„(x) = P exp \-i / gA-dy) (44) 




where the P-operator path-orders the AJs along the integration path. We 
use U^s in place of A^s on the lattice, because it is impossible to formulate 
a lattice version of QCD directly in terms of A^'s that has exact gauge 
invariance. The UJs, on the other hand, transform very simply under a 
gauge transformation: 

Upix) -> n(x)U fl (x)n(x + afi)K (45) 

This makes it easy to build a discrete theory with exact gauge invariance. 
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A link variable U„(x) is represented pictorially by a directed line from x 
to x + p,, where this line is the integration path for the line integral in the 
exponent of U^x): 



• • • 

x 



In the conjugate matrix U^x) the direction of the line integral is flipped 
and so we represent UUx) by a line going backwards from x + fi to x: 



• • • 

x 



A Wilson loop function, 

W(C) = ±Tr Pe- 1 ^^, (46) 

for any closed path C built of links on the lattice can be computed from the 
path-ordered product of the U^s and [H's associated with each link. For 
example, if C is the loop 



v 
t 



X 



then 



W(C) = i Tr (u^{x)U v {x + ap)... Ul(x) 



(47) 



Such quantities are obviously invariant under arbitrary gauge transforma- 
tions ©. 

You might wonder why we go to so much trouble to preserve gauge 
invariance when we quite willing give up Lorentz invariance, rotation in- 
variance, etc. The reason is quite practical. With gauge invariance, the 
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quark-gluon, three-gluon, and four-gluon couplings in QCD are all equal, 
and the bare gluon mass is zero. Without gauge invariance, each of these 
couplings must be tuned independently and a gluon mass introduced if one 
is to recover QCD. Tuning this many parameters in a numerical simulation is 
very expensive. This is not much of a problem in the classical theory, where 
approximate gauge invariance keeps the couplings approximately equal; but 
it is serious in the quantum theory because quantum fluctuations (loop- 
effects) renormalize the various couplings differently in the absence of exact 
gauge invariance. So while it is quite possible to formulate lattice QCD 
directly in terms of A^s, the resulting theory would have only approximate 
gauge invariance, and thus would be prohibitively expensive to simulate. 
Symmetries like Lorentz invariance can be given up with little cost because 
the symmetries of the lattice, though far less restrictive, are still sufficient 
to prevent the introduction of new interactions with new couplings (at least 
to lowest order in a). 

We must now build a lattice lagrangian from the link operators. We 
require that the lagrangian be gauge invariant, local, and symmetric with 
respect to axis interchanges (which is all that is left of Lorentz invariance). 
The most local nontrivial gauge invariant object one can build from the link 
operators is the "plaquette operator," which involves the product of link 
variables around the smallest square at site x in the fiu plane: 

Pfwix) = \ Re Tr (y^x)U u (x + a/i)E/J(x + aji + av)Ut{x)) . (48) 

To see what this object is, consider evaluating the plaquette centered about 
a point xq for a very smooth weak classical field. In this limit, 

iV « 1 (49) 

since 

Up « e~ l9aA » « 1. (50) 

Given that is slowly varying, its value anywhere on the plaquette should 
be accurately specified by its value and derivatives at x . Thus the cor- 
rections to (49) should be a polynomial in a with coefficients formed from 
gauge-invariant combinations of A^xq) and its derivatives: that is, 

Pyv = 1 ~ ci a 4 Tr (^^(xo)) 2 

- c 2 a 6 Tr (gF^(x Q )(Vl + B 2 u )gF^{x )) 

+ 0(a 8 ) (51) 
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where c\ and C2 are constants, and is the gauge-covariant derivative. The 
leading correction is order a 4 because F 2 is the lowest-dimension gauge- 
invariant combination of derivatives of A^, and it has dimension 4. (There 
are no F 3 terms because P^ v is invariant under — > or, equivalently, 
F -» -F.) 

It is straightforward to find the coefficients c\ and c 2 . We need only 
examine terms in the expansion of that are quadratic in A^; the cubic 
and quartic parts of F^ then follow automatically, by gauge invariance. 
Because of the trace, the path ordering is irrelevant to this order. Thus 



iReTrPe" 1 ^^ 



I Re Tr 



1 - i j gA ■ dx - \ (J> gA ■ dx*j + 0{A 3 ) 



(52) 



where, by Stoke's Theorem 

ra/2 



d>A-dx = / dx l _ l dx l/ [d l _ l A l/ (x + x) - d^A^xo + x)} 

Jn J-a/2 
ra/2 

= / dx^dxv [F^ v {xq) + ( x uPu + x^D fl )F llI/ (xo) + 

J-a/2 



a 2 F^(x ) + - (DJ + Bl)F^(x ) + 0(a , A 2 ). 



(53) 



Thus ci = 1/6 and c 2 = 1/72 in (H). 

The expansion in ( |5l|) is the classical analogue of an operator product ex- 
pansion. Using this expansion, we find that the traditional "Wilson action" 
for gluons on a lattice, 



(54) 



where (3 = 6/g 2 , has the correct limit for small lattice spacing up to correc- 
tions of order a 2 : 



S wn = J d 4 x^|iTrF^ + ^TrF^(D2+D^ + -.j. 



(55) 



We can cancel the a 2 error in the Wilson action by adding other Wilson 
loops. For example, the 2a x a "rectangle operator" 



R 



±Re Tr 



t 



-/' 



(56) 
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has expansion 

4 
6 



= 1 - - a 4 Tr (gF^) 2 ~^a 6 Tr (gi^(4Dj + D 2 v )gF t 



(57) 



The mix of a terms and a 6 terms in the rectangle is different from that in 
the plaquette. Therefore we can combine the two operators to obtain an 
improved classical lattice action that is accurate up to C(a 4 ) 13]: 



S, 



classical 



X,fl>U 



5P, 



[IV 



Rp,v ~\~ Run 

12 



+ const 



/d 4 x^iTrF^ + 0(c 



(58) 
(59) 



This process is the analogue of improving the derivatives in discretizations 
of non- gauge theories |] 

Exercise: The euclidean Green's function or propagator G for a nonrelativistic 
quark in a background gauge field satisfies the equation 



D 2 \ 

D t 

2M J 



G{x) = 5 4 (x) 



(60) 



where D p = <9 M — igA^(x) is the gauge-covariant derivative. Show that the 
static (M = oo) quark propagator is 



Goo(x,t) 



Pe 



-I f* gAo(x,t)dt 



5 3 (X) 



which on the lattice becomes 

Gocfat) = J7 t t (x,f-a)[/ t + (x,t-2a) ... C//(x,0). 



(61) 



(62) 



Propagation of a static antiquark is described by GJ^ and therefore the "static 
potential" V(r), which is the energy of a static quark and antiquark a dis- 
tance r apart, is obtained from 



W{r,t) = (0||Tr 



10} 



(63) 



■t- 



6 An important step that I have not discussed is to show that the gluon action is positive 
for any configuration of link variables. This guarantees that the classical ground state of 
the lattice action corresponds to F M „ = 0. See Q for a detailed discussion. 
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where for large t 



W(r, t) — > const e 



-V(r)t 



(64) 



Exercise: Defining the "twisted-rectangle operator' 




Re Tr 



(65) 



show that 



'classical 




+ const 



(66) 



d 4 z£±Tr^ + 0(a 4 ) 



(67) 



This is an alternative to the improved gluon action derived in the previous 



3.2 Quantum GluonsQ 

In the previous section we derived improved classical actions for gluons that 
are accurate through order a 4 . We now turn these into quantum actions. 
The most important step is to tadpole improve the action by dividing each 
link operator by the mean link no: for example, the action built of 
plaquette and rectangle operators becomes 



The no's cancel lattice tadpole contributions that otherwise would spoil 
weak-coupling perturbation theory in the lattice theory and undermine our 
procedure for improving the lattice discretization. Note that no ~ 3/4 when 
a = .4fm, and therefore the relative importance of the R^s is larger by 
a factor of 1/uq ~ 2 than without tadpole improvement. Without tadpole 
improvement, we cancel only half of the a 2 error. 

The mean link uq is computed numerically by guessing a value for use in 
the action, measuring the mean link in a simulation, and then readjusting 
the value used in the action accordingly. This tuning cycle converges rapidly 
to selfconsistent values, and can be done very quickly using small lattice 
volumes. The no's depend only on lattice spacing, and become equal to one 
as the lattice spacing vanishes. 

7 This section and the next are based on work with M. Alford, W. Dimm, G. Hockney 
and P. Mackenzie that is described in 1114] . 



exercise. 




(68) 
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The expectation value of the link operator is gauge dependent. Thus 
to minimize gauge artifacts, uq is commonly defined as the Landau-gauge 
expectation value, (0| | Tr ?7^|0)lg- Landau gauge is the axis-symmetric 
gauge that maximizes no, thereby minimizing the tadpole contribution; any 
tadpole contribution that is left in Landau gauge cannot be a gauge artifact. 
An alternative procedure is to define u$ as the fourth root of the plaquette 
expectation value, 

n = (0|P^|0) 1 / 4 . (69) 

This definition gives almost identical results and is more convenient for 
numerical work since gauge fixing is unnecessary. 

Tadpole improvement is the first step in a systematic procedure for im- 
proving the action. The next step is to add in renormalizations due to 
contributions from k>ix /a physics not already included in the tadpole im- 
provement. These renormalizations induce a 2 a s {ir /a) corrections, 



5C = asna 2 ^ Tr(F^BlF^) 
+ a s r 3 a 

Tr(D M F M<T D,F, CT ) 

+ (70) 

that must be removed. The last term is harmless; its coefficient can be set 
to zero by a change of field variable (in the path integral) of the form 

A^^ A^ + a 2 a s f(a s )J2VuF^. (71) 

Since changing integration variables does not change the value of an inte- 
gral, such field transformations must leave the physics unchanged.^] Opera- 
tors that can be removed by a field transformation are called "redundant." 
The other corrections are removed by renormalizing the coefficient of the 
rectangle operator R^ u in the action, and by adding an additional operator. 
One choice for the extra operator is 



C^va = I Re Tr 



(72) 

8 One must, of course, include the jacobian for the transformation in the transformed 
path integral. This contributes only in one-loop order and higher; it has no effect on 
tree-level calculations. 
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Then the action, correct up to 0(a 2 a 2 , a 4 ), is [15] 




} 




(73) 



where 



g 



g 



l + .48a s (vr/a) 
.055 a s (7r/a). 



(74) 
(75) 



The coefficients r g and c g are computed by "matching" physical quan- 
tities, like low-energy scattering amplitudes, computed using perturbation 
theory in the lattice theory with the analogous quantity in the continuum 
theory. The lattice result depends upon r g and c g ; these parameters are 
tuned until the lattice amplitude agrees with the continuum amplitude to 
the order in a and a s required: 



Note that tadpole improvement has a big effect on these coefficients. With- 
out tadpole improvement, r g = l+2a s ; that is, the coefficient of the radiative 
correction is four times larger. Tadpole improvement automatically sup- 
plies 75% of the one- loop contribution needed without improvement. Since 
a s ~ 0.3, the unimproved expansion for r g is not particularly convergent. 
However, with tadpole improvement, the one- loop correction is only about 
10-20% of r g . Indeed, for most current applications, one-loop corrections to 
tadpole-improved actions are negligible. 

Finally note that the twisted-rectangle action (|66|) must also be tadpole 
improved: 



The correction term in this action has two more uq's than that in our other 
action. This makes St r t much more dependent upon tadpole improvement. 
Comparing results generated by our two actions provides a sensitive test of 
tadpole improvement. The a s corrections to this action have not yet been 
computed, but are probably unnecessary. 

Exercise: Show that the gauge that maximizes (0|| Tr J7 M |0) becomes Landau 
gauge (d ■ A = 0) in the a— >0 limit. 




(76) 




(77) 
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Exercise: Tadpole diagrams are only important because they are ultraviolet diver- 
gent; the loop integrations generate factors of 1/a 2 that cancel the explicit o's 
from the link operator U^. Consequently such contributions can be signif- 
icantly reduced if the effective ultraviolet cutoff for the gluons is lowered. 
Lowering the cutoff by a factor of two, for example, would reduce tadpole 
contributions by a factor of four, rendering them harmless. This reduction is 
easily accomplished by smearing, where, for example, 

(4t-\4 \ n 
i-VJ Mx) - (78) 

The low-momentum components of the smeared A M field are identical, up 
to 0(a 4 ) corrections, with the unsmeared field; but at high momentum the 
smearing factor acts as an ultraviolet cutoff 

4 n4 \ n / 4 4 \ ™ 

1 ^ l- e -^-) «e-^ 4 . (79) 



Tadpole improvement should be unnecessary when the smeared field is used 
in quark actions, etc., in place of the original field. Parameter e is chosen to 
give the appropriate cutoff, while n is chosen sufficiently large that smearing 
factor approximately exponentiates as indicated. Design and test an imple- 
mentation of this smearing scheme for the link variables U^. 

Exercise: Lattice </> 4 theory has a 2 corrections of the form 

SC = c ia 2 Y, ( d l^f + c 2 a 2 f + c 3 a 2 (d 2 (f>) 2 + c 4 a 2 fd 2 <t>. (80) 

A* 

Find a field transformation 

<f>^i = <t> + a 2 f(<t>,%<f>...) (81) 

that removes the C3 and C4 terms (that is, the redundant operators). 
Note that <j> — > (j) + S(f> implies 

S[4>]^S + 54>^. (82) 

Thus the leading redundant operators can be rearranged so that they are 
proportional to 6S/S(f>. Since the classical field equation is SS/S(f> — 0, re- 
dundant operators are often said to have no effect "because of the equations 
of motion." More precisely it is because they can be removed by a change of 
integration variable in the path integral. 

Show that the jacobian \5(f)/5<p\ has no effect at tree-level, but must be in- 
cluded when quantum loop corrections are important. (Hint: Try computing 
a simple scattering amplitude using perturbation theory in the original and 
in the transformed theory.) 



25 



3.3 Tests of the Gluon Action 



Having a procedure for systematically improving the gluon action, we need 
to establish whether the improvements really do improve simulations. 

The effects of the 0(a 2 ) improvements in the gluon action are imme- 
diately apparent if one compares the static-quark potential computed, us- 



ing (|64|) , with the Wilson action (|54|) and with the improved action (173) [ 14 ] . 



Monte Carlo simulation results for these potentials are plotted in Fig. to- 
gether with the continuum potential. The Wilson action has errors as large 
as 40%. These are reduced to only 4-5% when the improved action is used 
instead. 

The dominant errors in the uncorrected simulation of V(r) reflect a fail- 
ure of rotational invariance, which is expected since the a 2 error in the 



Wilson action (54) is neither Lorentz nor rotationally invariant. The points 
at r = a, 2a, 3a . . . are for separations that are parallel to an axis of the 
lattice, while the points at r = \/2a, \^3a . . . are for diagonal separations be- 
tween the static quark and antiquark. In Table [l], I tabulate the error AV 
in the potential at r = \/3a for a variety gluon actions, with and without 
tadpole improvement and one-loop radiative corrections. As expected, the 
correction term in the action is significantly underestimated without tad- 
pole improvement. The tadpole-improved action is very accurate both with 
and without one-loop corrections, suggesting that (D(a 2 a s ) corrections are 
comparable to those of 0(a 4 ) and therefore unimportant for most current 
applications. 

I also include in Table [l] results obtained using the twisted-rectangle 
action St v t m (k7D- This action gives a potential that is essentially identical 
to that obtained with the other improved action. The a 4 errors are very 
different in the two actions; their agreement suggests that a errors are small. 
It also is striking confirmation of the importance of tadpole improvement 
since the uo's more than double the size of the correction term in the twisted- 
rectangle action when a = 0.4 fm. 

Another test of our improved action is to compute physical quantities on 
lattices with several different lattice spacings, looking for evidence of finite-a 
errors. These are consistently small. For example, the temperature at which 
gluonic QCD has a deconfining phase transition can be computed using the 
improved the action. One finds |l6[ |: 

r 298 (6) MeV for a = .32 fm 
T c = < 309(6) (6) MeV for a = .21 fm. (83) 
1 303 (8) MeV fora = .16fm. 
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Figure 4: Static-quark potential computed on 6 4 lattices with a ~ 0.4 fm 
using the Wilson action and the improved action. The dotted line is the 
standard infrared parameterization for the continuum potential, V(r) = 
Kr — 7r/12r + c, adjusted to fit the on-axis values of the potential. 
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Table 1: Error in the static quark potential at V(y/3a) for a variety of gluon 
actions. The lattice spacing in each case is a ~ Aim; K is the slope of the 
linear part of the static potential. 



Action AV(V3a)/KV3a 

unimproved (Wilson) .41 (2) 

tree-level improved, no tadpole improvement .15 (1) 

one- loop improved, no tadpole improvement .12 (2) 

tree-level improved, with tadpole improvement .05 (1) 

one-loop improved, with tadpole improvement .04 (1) 

twisted-rectangle correction, with tadpole improvement .04 (2) 



The results all agree to within statistical errors of a few percent. 

3.4 Anisotropic Lattices and Tuned Gluon Actions^ 

As a final illustration of improved gluon actions with tadpole-improved oper- 
ators, I now show results of simulations employing anisotropic lattices with 
temporal lattice spacings at that are smaller than the spatial lattice spac- 
ings a s . Such anisotropic lattices lead to greatly improved signal-to-noise 
in Monte Carlo simulations of hadrons; they are also useful for designing 
improved quark actions, as I discuss later. On an anisotropic lattice, our 
rectangle-improved action becomes 

c ^ a t (5P SS > 1 R ss , 1 R s > s \ 

x,s>s' 4 * 

where s and s' are spatial directions. The mean links 114 and u s are different 
on anisotropic lattices. When at < a s /2, u% is very close to its continuum 

9 This section is based on work with M. Alford, T. Klassen, C. Morningstar, M. Peardon 
and H. Trottier. 
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Figure 5: The static-quark potential computed on an anisotropic lattice in 
different orientations. Results are shown for at/a s equal 1/2. Except as 
shown, analysis errors are smaller than the plot symbols. The fit in each 
plot is to the open circles; the fitting function is V{r) = Kr — b/r + c. 

value of one. Thus we take 

u t = 1 (85) 
u s = (Pss>) 1/4 . (86) 

The same u s , to within a few percent, can be obtained from the link ex- 
pectation value in Coulomb gauge, which is the natural gauge choice for 
anisotropic lattices with small aj's. 

The anisotropic simulation can be tested by comparing the static quark 
potential computed from Wilson loops in different orientations: loops in the 
x—t plane, loops in the x-y plane, and loops in the t-x plane (with x playing 
the role of time). If the anisotropic theory is tadpole-improved all three 
potentials agree (Fig. Without tadpole- improvement the slopes disagree 
by 20% (Fig. §. 

The coefficients of the correction terms in the gluon action for anisotropic 
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Figure 6: The static-quark potential computed on an anisotropic lattice in 
different orientations, as in the previous figure. The action used here was 
not tadpole improved. The input value for a t /a s is 1/2. 
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lattices are known only to leading order in a s . While it is clear from the simu- 
lations that tadpole improvement captures the bulk of the renormalizations, 
these same simulations provide a handle on further radiative corrections. 
This is because the leading corrections in the gluon action affect the extent 
to which continuum symmetries, like rotation invariance, are restored. We 
can use symmetry restoration as a criterion for tuning the couplings beyond 
their tadpole-improved tree-level values. Thus, for example, the static po- 
tential at r = (a s , 2a s , 2a s ) should agree exactly with that at r = (3a s , 0, 0). 
One might tune the value of the spatial mean link u s used in the action 
until this is the case. Similarly one would like the anistropy at/a s built 
into the action to agree with the renormalized anisotropy coming out of the 
simulations. The ratio ut/u s can easily be tuned to make the two agree 
exactly. Such tuning does not capture all of the 0(a s ) corrections except 
in the simplest of situations. However it may be a useful intermediate step 
between tree-level tadpole improvement and a full perturbative analysis of 
the radiative corrections. 

Exercise: One defect of our improved gluon action is that the gluon spectrum has 
ghosts. These are caused by the R ts terms in the action, which extend two 
steps in the time direction. (Numerical ghosts, like these, are discussed in 
greater detail in the section on light quarks.) These gluon ghosts, being at 
high energies (of order 2.6/a t ), should have little effect on simulation results, 
particularly for smaller a t 's; this is confirmed by simulations. Nevertheless 
one might wish to remove them completely. This is easily done, when at is 
small compared with a s , by dropping the R ts terms in the action. This 
introduces new errors of order a\ but these will be negligible relative to the 
a\ errors provided at is small enough. Modify the st terms in the action so 
that all a 2 s errors are still canceled in the absence of the R ts contribution. 

3.5 Summary 

In this section I have outlined how one designs actions for simulating gluon 
dynamics on coarse lattices. The simulation results show that accurate 
simulations are possible even on lattices as coarse as a = .4 fm using very 
simple actions. Tadpole improvement is essential, but one-loop radiative 
corrections to the action are not too important. Rotational invariance of 
the potential and, for anisotropic lattices, space-time interchange symmetry 
provide sensitive tests of improved gluon actions, and can be used to tune 
the leading correction terms. 

The coarse lattices make simulations very fast. Most of the simulation 
results shown here were obtained using an IBM RS6000/250 desktop work- 
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station, which is powered by a personal-computer CPU (66MHz PowerPC). 
Any of these results could be generated easily on a high-end personal com- 
puter. 



4 Light Quarks on Coarse LatticesQ 

In this section I develop techniques for designing very accurate discretiza- 
tions of the Dirac equation for quarks in a gluonic field. I begin by reviewing 
various resolutions of the "doubling problem" which plagues even the sim- 
plest discretization. I then discuss two currently popular discretizations, 
and finally new discretizations that are significantly more accurate. 



4.1 Naive Discretization and Doubling 

The euclidean Dirac lagrangian in the continuum is 

£ = ^(D-7 + m)V (87) 



where 



and I take 



1 \ a s \ _ . / 1 \ 

^={0 -1 J ls= [a s J 75 =-WW& = * [ _x oj 

(89) 

where, again, s and t as subscripts signify spatial and temporal indices 
respectively. The obvious discretization is 

Aat = V> (A • 7 + m + 0(a 2 )) i/i (90) 

where now A^ is a (tadpole-improved) gauge-covariant derivative: 

a / / \ U pk {x)^(x + a[i)-Ul{x-a[i)' l lj{x-aji) 

&H^{x) = 7T- ~ • (91) 

lauQ 

The other finite-difference operators, A^ , A^ . . . , can all be made gauge- 
covariant by adding U^s as in A^; in what follows I will always use these 
symbols to refer to the gauge-covariant differences. 

10 This section is based on work with M. Alford and T. Klassen. 
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Unfortunately this simple lattice Dirac action has a very serious pathol- 
ogy To see this consider the p = spin-up solution of the corresponding 
Dirac equation without gluon fields: 



(92) 



(There is no I in the exponent because t is euclidean; E is the minkowski 
energy) Substituting into the Dirac equation, 

(A-7 + m) Vt(*) = 0, (93) 

we obtain 

e —aE p aE 

(A t + m) e~ Et = =► + m = 0. (94) 

Z. Ob 

This equation is quadratic in z = exp(aE) and therefore it has two solutions 
for E. We only want one! The two solutions are 

_ ( m + O(o?) (normal) , , 

\-m + Iir/a + 0(a 2 ) (doubler) " [ ' 

The first is the solution we want; the second is known as a "doubler." 

The possibility of a second solution for the energy is obvious from in- 
spection of the lattice Dirac equation. The problem lies in the A t 4> term. 
This links ip's two lattice spacings apart which means that the dispersion 
relation is second order in z = exp(aE). In general one expects n roots if 
the V>'s are spread over n lattice spacings — a potential problem when we go 
to improve the discretization. The doubler cannot be ignored. It is related 
to the normal solution by a symmetry of the gluon-free lattice theory: if 
V>(x, t) is a solution of the lattice Dirac equation then so is 

^(x,0^7 57t (-l) t/a V(x,t). (96) 

This symmetry operation turns the normal solution into the doubler solu- 
tion. Consequently the doubler acts as an extra flavor of quark, which is 
degenerate with the normal quark (when there are no gauge fields) . Further- 
more, in euclidean space there is nothing unique about the time direction. 



-Et 
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Thus there is an analogous doubling symmetry for each direction, and suc- 
cessive transformations involving two or more directions also are symmetries. 
This implies that there are fifteen doublers in all. 

The root of the doubler problem lies in the finite-difference approxima- 
tion for first-order derivatives. Fourier transforming d and A we obtain 

d -> lp (97) 
A -► lsm(pa)/a. (98) 

These agree well for low momenta, but, as p approaches the ultraviolet 
cutoff ir /a, the transform of A vanishes rather than becoming large. As a 
consequence, for example, high-momentum states can have very low or even 
vanishing energies. This is disastrous. Note that the same problem does not 
arise for second-order derivatives: 

8 2 -p 2 (99) 

A (2) - (- sin(pa/2)^ . (100) 

Although the transform of A^ 2 ) becomes quite different from p 2 at large p, 
it nevertheless remains large, going to its maximum of 4/a 2 at the ultravi- 
olet cutoff. Second-order (and higher even-order) derivatives do not cause 
doubling problems. 

I will discuss two approaches to the doubling problem. One involves 
adding operators to the quark action that are almost redundant but remove 
the doublers or drive them to high energies. The other involves "staggering" 
the quark degrees of freedom on the lattice. 

Exercise: Show that gauge transformation 

tp(x) -» fi(ar) i/>(x) U^x) -» 9,{x)U ll {x) Q(x + afi)' ! . (101) 

implies that 

A^(aO->ft(aOA„^(aO. (102) 
Also show that (setting Uo = l) 

p aD M _ p-aD^ 

A, = 2a ■ (103) 

(Hint: Why can't there be F^'s on the right-hand side?) This relation 
is crucial to whot follows; find a different definition of the the covariant 
difference operator A for which this relation is incorrect in 0(a). 
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4.2 Field Transformations to Remove Doublers 



The first approach to removing the doublers is to introduce into the action an 
operator that is redundant up to 0{a 2 ) and has higher order non-redundant 
parts that drive the doubler to high energies. A redundant operator remem- 
ber, is one that can be eliminated from the action by a change of variables 
on the quark fields in the path integral, and hence has no effect on any 
physical observables. For example, consider the field transformation from 
the continuum field ip c to a new field ip where 



1 



r a 



(A • 7 - m c ) 



and 



With this transformation our naive lattice action becomes 
ij) c (A • 7 + m c ) ip c = ip Q (A • 7 + m c ) if) 



(104) 
(105) 

(106) 



r a l . , a ■ gF 



A-7 + m- — A-A + 



i) + 0(a 2 ) (107) 



where m = m c + ram 2 /2. In deriving the final formula we used the lattice 
version of the identity 

a-gF 



(D- 7 ) z = D-D + 



where 



a 



i:.U — 2 [TA" T"] 



(108) 



(109) 



The new 0{a) term in the transformed lagrangian does not introduce new 
0{a) errors because it results from a field transformation; it is redundant 
in the sense discussed above for gluon operators. In its present form this 
term also has no effect on the doubling problem. The trick is to replace the 
operator A • A, which vanishes at the ultraviolet cutoff, with A^ 2 ), which 
is the same up to 0{a 2 ) errors in the infrared but large at the ultraviolet 
cutoff. The A^ 2 ) pushes the doublers off to high energies, of order 1/a, where 
they are harmless. Thus an C(a)-accurate quark action that is doubler free 
is 



£sw = V> 



A • 7 + m — 



r a 



a-gF 



2 V 2 
^ c (D- 7 + m c )Vc + 0(a 2 ), 



(110) 
(111) 
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where, again, 

i(j c = Qip m = m c + ram 2 /2. (H2) 

This is the Sheikholeslami-Wohlert (SW) quark action IT?]]. It becomes Wil- 
son's original action if the a ■ gF term is dropped (resulting in 0(a 1 ) errors). 

We can check that the doublers are removed by again substituting the 
p = solution into the corresponding Dirac equation. We obtain 

A t + m -^A( 2 ))e-^ = 

(1 - r) e~ aE - (1 + r) e aE + 2r + 2 am = 0. (113) 

which has roots 

F = ( m c + 0(a 2 ) (normal) 
\(l/a)ln[(r-l)/(r + l)] + -" (s host ) ' 

The doubler has been replaced by a high-energy "ghost" particle, as desired. 
The ghost goes away completely for r = l. 

Exercise: Create a lattice version of F M „ for use in the SW quark action. Build 
your operator out of tadpole-improved link operators. Try designing a dis- 
cretization that is accurate up to corrections of 0(a i ). 

4.3 a 3 -Accurate Lattice Dirac Theory 

Our numerical experiments in the first Section suggest that the a 2 errors in 
the SW quark action might still be quite large when the lattice spacing is of 
order .4 fm. Thus it is worthwhile developing a quark action that has errors 
only of order a 4 and higher. 

A attempt at such an improved Dirac theory is the action 



a 2 



A • 7 — — A (3) • 7 + m t 
6 



i/> (115) 



where the a 2 correction cancels the leading errors induced by the A • 7 term. 
Not surprisingly this action is plagued with doublers and ghosts. Since the 
temporal derivatives extend two steps in each direction, the energy disper- 
sion relation has four branches instead of one. 

There are several approaches one can take to solving this doubler disas- 



ter. The "D234" quark actions rely upon three ingredients 18]: 
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1. anisotropic lattices, with at<a s /2, which push ghost states to higher 
energies, where they decouple. This works because ghost energies are 
governed by the temporal lattice spacing and grow like 1 /at when at is 
reduced. Reducing the temporal lattice spacing also makes Monte 
Carlo measurements of hadron propagators much more efficient. 

2. redundant operators like 



in the quark action, which can be tuned to remove some of the doublers 
and ghosts. 

3. "irrelevant" operators like 



in the quark action, which have negligible effect on low-momentum 
physics (if a t < 02/2) but are useful at high energies for removing 
doublers and ghosts. Terms involving higher-order spatial derivatives, 
like a^A^\ can also be added. 

A variety of actions can be created by combining differing amounts of redun- 
dant and irrelevant operators. In general one tries to tune these additions 
to minimize the number and importance of ghost states, to give the best 
low-momentum behavior, and to minimize pathologies at high momentum. 

A simple example of a D234-type action is obtained by generalizing our 
derivation of the SW action. We begin with a field transformation ip c — ► ip 
where 




(116) 




(117) 



■ip c = Slip 



(118) 



with 




(119) 



(120) 



and 



(121) 




(122) 
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Applying this transformation to the naive improved action ( |115| ) and adding 
an irrelevant operator, we obtain^ 



£D234x 



m 



ra t 
2 



(123) 



This action is the same as the continuum action up to corrections of order a% 
and cig . It is easy to show that the choice 



1 r 
12 ~ 24 



(124) 



removes one of the ghost /doubler states, leaving two ghosts that have high 
energies for any value of r near 1. If in addition r = 2/3 is chosen, a second 
ghost disappears and only one remains (with an energy of order 2 /at-) The 
free-quark dispersion relation for this set of parameters is shown in Fig. |7] 
for a lattice with a s = 2 at = .4 fm; the ghost is roughly 2 GeV above the real 
quark and so harmless in most applications. 



4.4 Tests of Quark Action 

Highly improved actions like the D234x action ( |123[) have a large number of 
bare coupling constants that must be computed in perturbation theory — 
one for each term in the action. In the previous section we worked only to 
lowest order (tree-level) in perturbation theory. However our discussion of 
the gluon action suggests that this is sufficiently accurate provided all opera- 
tors are tadpole improved; the renormalizations due to k>ir/a physics affect 
phyical quantities by only a few percent. How can we ascertain whether the 
same is true of our quark actions? 

11 The designers of this action, having examined a large number of similar actions, all 
very effective, are currently grappling with a nomenclature crisis. I refer to the particular 
action shown here as the "D234x" action, but this is certainly not its final name. Most of 
what I say about this action is true of almost all other D234 actions, and so "D234x" can 
usually be taken to refer to a generic D234 action. 
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Figure 7: Dispersion relation for a free massless lattice quarks described by 
the D234x action. A low-energy, physical branch and a high-energy ghost 
branch are shown, together with the continuum dispersion relation (dotted 
line). Energies are for momenta proportional to (1, 1,0). 
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One powerful approach is to verify that continuum symmetries are re- 
stored by the corrections. For example, the leading error in the discretized 
Dirac operator, A-7 + m, breaks Lorentz invariance. The a^A^ 7^ term in 
the D234x action is meant to cancel this error. We can check that this term 
has done its job by comparing the dispersion relation of a meson (or baryon) 
computed in a simulation with the Lorentz-invariant dispersion relation of 
the continuum. In a Lorentz-invariant theory, a hadron's energy E^(p) and 
its three momentum p are such that the combination 

c 2 (p) = < (125) 

is independent of p and equal to the square of the speed of light (= 1). In 
Fig. ^ I show D234x simulation results for meson's C(p). The lattice for this 
simulation is very coarse, with a s = 2 at = .4 fm, and yet c(p) is within a 
few percent of one all the way out to momenta of order 1 GeV (~ 2/a s ). 
For comparison I include the same quantity calculated using the SW action 
(on an isotropic lattice), which has no correction for the o 2 error; it has 
25% errors by p= 1 GeV, and the momentum dependence of the discrepancy 
is consistent with the quadratic errors expected. This comparison shows 
that the coefficient of the of the a 2 correction in the D234x action is not 
renormalized by more than 10-20%, and that this renormalization will have 
no more than a few percent effect on low-momentum physics. The momenta 
used for the Monte Carlo points shown in this plot include momenta that 
are parallel to a lattice axis, and momenta that point along lattice diagonals; 
the fact that different momenta agree to within a few percent shows that 
rotation invariance is also restored by the correction term. 

The c(p) plot provides a second important check on the D234x action. 
On anisotropic lattices, the symmetry between space and time is not exact. 
Consequently the spatial and temporal parts of the leading Dirac operator A- 
7 are renormalized differently: 

D-7^ A t7t + c A-7, (126) 

where cq is a bare speed-of-light, 

Co = 1 + const a s + ■ ■ ■ . (127) 

Such a renormalization would shift the low-momentum values of c(p) away 
from one. In fact, with only our tree-level action, c(p) = 1.025(20) at 
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c(p) 0.6 - 




Figure 8: The speed of light computed using the D234x and SW quark 
actions from the energy and three momenta of a meson. The fit to the SW 
data is cq + k (ap) 2 where cq and k are constants. 
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Table 2: Mass differences between two simulations, one with r = 1 and the 
other with r = 2/3, for different quark masses. The bare quark masses in 
the two simulations were tuned to give the same m n /m p . 







Am p 






.65 


0.6 (2) % 


-0.13(13)% 


-0.3 (2) % 


0.4 (3) % 


.57 


0.6 (4) % 


-0.6(1.0)% 







low momenta; parameter cq must be within a percent of one. Note that 
this success depends crucially on tadpole improvement; without tadpole 
improvement c(p) = 1.2 at low momenta on our a s = .4fm lattice. 

The only other potentially large renormalization is in the 0(rat) terms 
induced by the field transformation. Again there is a symmetry associated 
with this term: field-redefinition or "r symmetry" . The operator multiplied 
by r is redundant and so should have no effect on the low- momentum physics. 
Consequently the low-momentum physics should be invariant under changes 
in the value of r. If the renormalizations of the A^ 2 ) and a ■ gF terms in the 
action are very different, then the tree-level operator is no longer redundant 
and results must change as r is varied; on the other hand, if the renormal- 
izations are roughly the same then the net effect of the renormalizations is a 
harmless shift in r.P| The results of such an r-test for a close relative of the 
D234x action are shown in Table [2|. These show that the pion, rho, nucleon 
and delta masses all shift by less than a percent when r is changed from 2/3 
to 1. Since omitting the a- gF term shifts masses by 20%, our r-test implies 
that that the relative renormalization given by tree-level tadpole improve- 
ment is correct to within roughly 10-20%. A 10% change in the coefficient 
of a ■ gF would shift the meson masses by only 2%. 

In [19] a different test of the a ■ gF coefficient is proposed. This is based 
on the fact that the r terms in the action break the quark's chiral symmetry 
even for massless quarks. If the a ■ gF coefficient is properly tuned, the 
chiral symmetry breaking terms are redundant and the chiral symmetry 
Ward identity should be restored. Using Monte Carlo evaluations of the 
various terms involved, the authors tune the coefficient until the identity 

12 I am simplifying the discussion here. On an anisotropic lattice, where space-time 
interchange is not an exact symmetry, the temporal and spatial parts of A' 2 ' and a ■ gF 
have different radiative corrections. Thus the r-test is checking on the relative coefficients 
of four operators, as opposed to two in the isotropic case. 
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is satisfied. Working only with the SW action and relatively small lattice 
spacings, they find renormalizations that are larger than those above for 
D234: the coefficient of a ■ gF in the SW action is a little more 20% larger 
than the tadpole-improved tree-level coefficient at a~.lfm. Such a shift at 
such small lattice spacings has only a small effect on hadron masses; but in 
the same study the authors encountered pathological behavior in the SW 
action, with this correction, at large lattice spacings. Also the correction 
to tadpole improvement, though small throughout the range of their study, 
grows much more rapidly than naively expected as a increases from .03 fm 
to .1 fm. It would be very interesting to verify these results using the r-test, 
which provides an independent method of testing the same correction. 

Our tests of the D234 action confirm that tree-level improvement of the 
quark action, using tadpole-improved operators, is quite accurate. However 
the work I just described concerning the SW action suggests that nonpertur- 
bative tuning of parameters in the action may be important for some actions 
and operators. Such tuning is straightforward for D234: the same tests I 
described above can be used to tune the leading corrections in the action. In 
particular, the cq parameter, which is potentially the most important, can 
be easily adjusted so that c(p = 0) = 1. 

Finally, in Fig. ^, I show simulation results for the rho mass as a function 
of lattice spacing. I give results for the D234x action 1 15], which should 
have errors, for the SW action ||(|, which should have a 2 errors, and 
for the Wilson action [21], which should have a 1 errors. The D234x and 
SW simulations use an improved gluon action; the Wilson results are for 
an uncorrected gluon action. Assuming that the a dependence of the errors 
is as expected, all three actions agree quite well on the continuum limit 
(a — > 0) mass of the rho, but the more improved the quark action, the 
sooner it gets there. This figure demonstrates very clearly that we know 
how to systematically improve quark actions for use on very coarse lattices. 
Furthermore it is clear that simulations that are accurate to within a few 
percent are possible on lattices with spatial lattice spacings between .3 and 
.4fm. 

4.5 t-Staggered Quarks 

A second approach to removing doublers is to "stagger" the quark degrees 
of freedom on the lattice. If we rewrite the Dirac field in terms of two 
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Figure 9: Rho mass versus lattice spacing computed using the Wilson quark 
action (W), the Sheikholeslami-Wohlert action (SW) and the D234x action. 
The charmonium S-P splitting was used to determine a s . These simulations 
have no quark vacuum polarization and therefore the rho mass is expected 
to differ from its experimental value. 
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two-component fields, (ft and x> 



i> = (it) t/at y_ x (128) 

the naive lattice Dirac equation (|9(]) falls apart into two decoupled two- 
component equations: 

(m c + A t + (-l)*/ at A-<x)</> = (129) 

(m c + A t + (-l) 1+ */ a 'A-cr) x = (130) 

The second equation is identical to the first equation, but shifted one lattice 
spacing in the t direction; it leads to identical quark physics. One of these 
equations describes the quark while the other describes the degenerate dou- 
bler. We get rid of the (temporal) doubler simply by dropping one of the 
two equations, say the second one. The remaining equation describes two 
quark polarizations and two antiquark polarizations, but no doublers. This 
approach is called "staggered" because dropping the x field is equivalent to 
specifying only half the spinor components of the original Dirac field on time 
slices of the lattice with even t/at, while only the other half are specified 
on time slices with odd t/at. The spinor components are staggered between 
time slices. 

This discretization still has a 2 errors and the spatial derivatives have 
poor behavior at high momentum. These problems are easily remedied by 
replacing 

n 2 

A -» A-^A® (131) 

,5 



m — ► m — so, 



(A( 2 )) (132) 



in (129), where the correction terms involve spatial derivatives only. The 
parameter s is set at some small value (~ .02) that is large enough to keep 
the high-momentum energy large but small enough that is has little effect 
on the low- momentum physics. The corrected equations still have a 2 errors, 
but these are negligible on anisotropic lattices with at < a s /2. 

This action has never been tried. It should be as accurate as the D234- 
type actions. It has half as many fields, no field transformation, no ghosts, 
and is somewhat simpler than the D234 actions; it may well be competitive. 

A more conventional action, the Kogut-Susskind or staggered-quark ac- 
tion, is obtained by staggering the quark degrees of freedom in spatial as 
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well as temporal directions. This reduces the number of doublers to only 
three (from fifteen in the original naive action). Rather than removing the 
remaining doublers, in this approach one treats them as three extra flavors 
of quark. If these extra flavors were degenerate with the original quark, this 
approach might be computationally very efficient, perhaps giving four times 
the statistics for the same computational cost as other discretizations. Un- 
fortunately the gluonic interactions lift the degeneracy between the flavors 
in the simplest version of this theory, resulting in significant mass differ- 
ences. It would be very interesting to develop improved versions of this 
theory where the degeneracy is more accurately preserved. 



4.6 Summary 

Temporal derivatives usually cause problems when improving the discretiza- 
tion of an action. This is because successive improvements to the tempo- 
ral derivatives involve operators that are increasingly nonlocal in time and 
therefore that almost inevitably result in ghosts of one sort or another. This 
problem is unusually serious for quark actions where even the leading-order 
discretization leads to doublers. 

As illustrated in the previous sections, there are two standard techniques 
for removing the doublers. One is to introduce redundant operators that 
break the doubling symmetry; this leads to the SW and D234-type actions. 
The second is to stagger the quark spinor components on the lattice; this 
leads to actions for f-staggered quarks and Kogut-Susskind quarks. 

Having repaired the leading-order operator, we face an even more acute 
problem in the af correction to Atjt- This can be dealt with, as in the D234 
actions, by using anisotropic lattices, with reduced aj's, to push the ghosts 
to high energies, and by adding irrelevant operators to cancel some of the 
ghosts. Alternatively, using anisotropic lattices, one might leave the af er- 
rors uncorrected; these will be negligible compared with errors provided 
at is small enough compared with a s . (Numerical experiments suggest that 
the a\ errors are only a few percent when at is less than .2fm.) A more 
complicated alternative, which I have not discussed, is to add irrelevant 
operators that replace the a^/^f 1 ^ correction term with less troublesome 
operators: for example, one can replace 



Aj 3) 7t -\ [(A • 7 + m c ) (A t7t ) z + (A t7t )" (A • 7 + m c )j (133) 

-► (A • 7 + m c ) Aat (A • 7 + m c ) (134) 
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7 + m c , m c 



— a 



■9F/2} 



(135) 



using successive field transformations and obtain, in the last step, an oper- 
ator that has only spatial derivatives. 

The numerical results in this section suggest that these various strategies 
work about as well as expected, and that accurate simulations are possible 
with spatial lattice spacings between .3 and .4fm. They also show that the 
meson dispersion relation (i.e., c(p)) at high and low momenta, the r-test 
and the p mass are all very sensitive to finite-a errors. 

5 Heavy Quarks 

As a final illustration of improved actions I briefly review techniques for 
simulating heavy quarks like the c and b quarks. Heavy quarks play a cen- 
tral role in the experimental study of weak interactions, and consequently 
it is important that we have reliable simulation techniques for their study. 
However heavy quark dynamics is not easily simulated using the quark ac- 
tions of the previous section. Those actions only work well for quarks with 
energies and momenta small compared with the ultraviolet cutoff ir/a. But 
heavy-quark hadrons have large rest energies (masses) and so necessitate 
small lattice spacings. Simulating an T, for example, using such techniques 
requires a lattice spacing of order 1/9.4 GeV = .02 fm or less — much too 
costly with today's computers. Fortunately the heavy quarks are generally 
nonrelativistic in such mesons. Most of such a meson's energy is due to the 
rest mass energy of its heavy-quark constituents, and rest mass energy plays 
little role in the dynamics of nonrelativistic particles. Consequently we can 
develop accurate simulation techniques that work well when the lattice spac- 
ing is chosen according to the size, rather than the mass, of the heavy-quark 
hadron: that is, for example, a = .lfm rather than .02 fm for upsilonsj 1 ^ 

Here I discuss only the NRQCD approach to heavy quark dynamics, 
where one takes maximum advantage of the fact that the heavy quarks 
are nonrelativistic [ p3| ] . I illustrate the technique by reviewing simulation 
results for the upsilon and psi families of mesons. These are attractive 
systems to study because so much is known already about them. The quark 

13 Actually since it is only the energy that is large, only the temporal lattice spacing 
need be small. Consequently the actions from the previous section can be used efficiently 
if one works with anisotropic lattices and small a t 's. Experiments using D234 actions 
are very encouraging Jl8[ . Yet another approach to heavy-quark dynamics can be found 
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potential model, for example, provides an accurate phenomenological model 
of the internal structure of the mesons. Also there is much experimental 
data for these mesons. The low-lying states are largely insensitive to light- 
quark vacuum polarization — for example, the T, T', Xh ■ ■ ■ are all far below 
the threshold for decays into B mesons — and therefore can be accurately 
simulated with rif = light-quark flavors. Furthermore, these mesons are 
very small; for example, the T is about five times smaller than a light hadron. 
This makes them an excellent testing ground for our ideas concerning large 
lattice spacings. 

Given that the heavy quarks are nonrelavistic, the most important mo- 
mentum scales governing a heavy-quark meson's dynamics are smaller than 
the heavy quark's mass. We can take advantage of this fact by choosing an 
inverse lattice spacing of order the quark mass, thereby excluding relativis- 
ts states from the theory. Then it is efficient to analyze the heavy-quark 
dynamics using a nonrelativistic lagrangian (NRQCD). The lagrangian used 
to generate the results I show below was 

£nrqcd = E{^n^+^(l-^)V/(l-^) n (l-a t <5i7)^(x 

- -£tf(x)il>(x), (136) 

X 

where n = 2, Hq is the nonrelativistic kinetic-energy operator, 

A {2) 

H ° = < 137 » 

Mq is the bare heavy-quark mass, and 5H is the leading relativistic and 
finite-lattice-spacing correction, 

8M 3 V 2n J 24M 
1 ■<r-B + J* (A-E-E-A) 



2M 8M n 2 



o 

^■(AxE-ExA). (138) 



8M 



Here E and B are the chromoelectric and chromomagnetic fields. As usual, 
the entire action is tadpole improved by dividing every link operator by 
uq. Potential models indicate that corrections beyond SH contribute only 
of order 5 MeV to T energies, and two or three times this for ^'s. 
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The details of this lagrangian are unimportant to us here. What I want 
to focus on is the extent to which the improvement program outlined in 
earlier sections is effective for heavy-quark dynamics. The 5H term consists 
of corrections just like the ones we have been analyzing for the other parts of 
QCD, the only difference here being that we are correcting both for finite-a 
and for the absence of relativity (ie, order v 2 /c 2 errors). The coefficients 
of the correction terms were determined with tree-level perturbation theory 
and tadpole improvement, using precisely the techniques outlined in the 
earlier sections. The extent to which 5H improves the simulation results is 
a measure of the efficacy of all the techniques discussed in these lectures. 

Consider first a set of simulations of the T family that used a lattice spac- 
ing of about 1/12 fm — roughly half the meson's radius [24]. The spectrum is 
very well described by the simulations; results for the low-lying spectrum in 



different channels are shown in Figure 10. These compare well with experi 



mental results (the horizontal lines), as they should since systematic errors 
are estimated to be less than 10-20 MeV. Two sets of Monte Carlo results 
are shown: one has no light-quark vacuum polarization, the other includes 
u and d quarks; as expected quark vacuum polarization has little effect on 
the specturm. It is important to realize that these are calculations from first 
principles. The only inputs are the lagrangians describing gluon and quark 
dynamics, and the only parameters are the bare coupling constant and the 
bare quark mass. In particular, these results are not based on a phenomeno- 
logical quark-potential model. These are among the most accurate lattice 
results to date. 

The corrections in 5H have only a small effect on the overall spectrum, 
but the spin structure is strongly affected. The lagrangian without 5H is 
spin independent, and gives no spin splittings at all. Simulation results 
for the spin structure of the lowest lying P state are shown in Figure [U]. 
Again these compare very well with the data, giving strong evidence that 
corrected lagrangians work. Systematic errors here are estimated to be of 
order 5 MeV. Note that the spin terms in 5H all involve either chromoelectric 
or chromomagnetic field operators. These operators are built from products 
of four-link operators and so tadpole improvement increases their magnitude 
by almost a factor of two at the lattice spacing used here. Simulations 
without tadpole improvement give spin splittings that are much too small. 

The correction terms also affect the extent to which the simulation repro- 
duces continuum symmetries. For example, the low-momentum dispersion 
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Figure 10: NRQCD simulation results for the spectrum of the T( 3 Si) and 
an d their radial excitations. Experimental values (dashed lines) are 
indicated for the S-states, and for the spin-average of the P-states. Simula- 
tion results are for rif = (solid circles) and rif = 2 (open circles) light-quark 
flavors. The energy zero for the simulation results is adjusted to give the 
correct mass to the T. 
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Figure 11: Simulation results for the spin structure of the lowest lying P- 
state in the T family. The dashed lines are the experimental values for the 
triplet states. Simulation results are for rtf = (solid circles) and n f =2 (open 
circles) light-quark flavors. 



relation of an upsilon in the simulation must have the form 

£ T (p) = M 1 + ^- + ---. (139) 

In a relativistically invariant theory M\ =M<i is the upsilon mass. However 
in a purely nonrelativistic theory M<i equals the sum of the quark masses, 
which differs from M\ since M\ also includes the quarks' binding energy. 
Only when relativistic corrections are added is Mi shifted to include the 
binding energy. Simulations tuned to the correct mass M\ = 9.5(1) GeV 
give M2 = 8.2(1) GeV when 5H is omitted. When 5H is included, the 
simulations give M% = 9-5(1) GeV, in excellent agreement with M%. All of 
the spin-independent pieces of 8H contribute to the shift in M2; once again 
we have striking evidence that corrected actions work. 

This T simulation has only two parameters: the bare coupling constant 
and the bare quark mass. These were tuned until the simulation agreed with 
experimental data. From the bare parameters we can compute the renor- 
malized coupling and mass. This simulation implies that the renormalized 
or "pole mass" of the b quark is [25] 



M b = 5.0 (2) GeV. (140) 
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Figure 12: S, P and D states of charmonium computed on lattices with a 
range of lattice spacings. Improved quark and gluon actions were used in 
each case except a = .17fm where only the quark action was improved. 



The renormalized coupling that is obtained corresponds to [ 



«§(M Z ) 



.1175(25) from Xb~~f splitting 
.1180 (27) from T'-T splitting 



(141) 



depending upon which T splitting in the simulation is tuned to agree ex- 
actly with experiment. These values agree with results from high-energy 
phenomenology (but are more accurate). This last result is striking: it 
shows that the QCD of hadronic structure and the QCD of high-energy 
quark and gluon jets are really the same theory. 

Another check on the corrected action is to show that results for physical 



quantities are independent of the lattice spacing. In Fig. |12| I show the low- 
lying spin-averaged spectrum of the i/j mesons computed using very coarse 
lattices, with a ranging from .17 fm to .40fm|]lJ]. The IS and IP masses 
are tuned; the 2S and ID masses are predicted. The latter masses are 
independent of lattice spacing to within statistical errors of a few percent. 



In Fig. 13 I show the IS and IP radial wavefunctions for charmonium, 
computed on lattices with two different lattice spacings using the (improved) 
NRQCD actions [jl4| . The most striking observation from these pictures 
is that the charge radius of the ip is almost exactly equal to one lattice 
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spacing on the coarser lattice, and yet the results from the coarser lattice 
are essentially identical to those from the finer lattice. These data confirm 
that in general one needs only a few lattice points per direction within a 
hadron to obtain accurate results (few percent) from an a 2 -accurate action. 
Numerical experiments with the NRQCD action, where a 2 corrections are 
easily turned on and off, suggest that 

(0{a?) corrections ^ 15-20% 

Hiadron-a I (142) 

lo(a s a 2 ,a 4 ) ^2-4% 

NRQCD simulations of heavy quarks were among the very first in which 
tadpole improvement played a crucial role. They remain among the best 
demonstrations of the technique. 



Exercise: Verify that the continuum limit of the lattice NRQCD lagrangian (136) 
is a Schrodinger action with the standard relativistic corrections (chromo- 
magnetic moment, Darwin term, spin-orbit coupling. . . ). 



6 Conclusions 



Improved discretizations of QCD dynamics and large lattice spacings work. 
Tree- level (i.e., classical) improvement of the quark and gluon actions give 
results that are accurate to within a few to several percent for lattice spac- 
ings as large as .3-.4fm, provided the all operators are tadpole-improved. 
Tree-level improvement is not surprisingly effective; it isn't perfect. Rather, 
it works about as well as one would naively expect: there are certainly 0{a s ) 
corrections (and probably also nonperturbative corrections), and these will 
shift couplings by 10 or 20% and physical quantities by a few percent. That 
QCD simulations now conform to naive expectations and intuition is a ma- 
jor advance, perhaps the major advance. The current level of precision is 
more than sufficient to have a large impact on QCD phenomenology. And 
these techniques provide a firm foundation, a starting point, for future high- 
precision work that will rely upon some mixture of systematic perturbative 
calculations of the corrections, nonperturbative tuning, nontrivial transfor- 
mations of the fields like those that lead to "perfect actions", and perhaps 
analytic techniques like strong-coupling expansions. 

The potential impact on QCD simulations of coarse lattices is enormous. 
This potential will be realized if lots of people become involved in the design, 
implementation and application this technology. 
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Figure 13: The radial wavefunctions for the IS and IP charmonium com- 
puted using improved actions and two different lattice spacings. Wavefunc- 
tions from a continuum quark model are also shown. Statistical errors are 
negligible for the 15 wavefunction. 
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